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ABSTRACT 

>' 

t^j- t We present a preliminary analysis of the X-ray spectrum of the accretion disk corona 

source, 4U 1822—37, obtained with the High Energy Transmission Grating Spectrom- 
eter onboard the Chandra X-ray Observatory. We detect discrete emission lines from 
photoionized iron, silicon, magnesium, neon, and oxygen, as well as a bright iron flu- 
orescence line. Phase-resolved spectroscopy suggests that the recombination emission 
comes from an X-ray illuminated bulge located at the predicted point of impact between 
the disk and the accretion stream. The fluorescence emission originates in an extended 

region on the disk that is illuminated by light scattered from the corona. 
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Subject headings: stars: individual (4U 1822—37) — stars: binaries: eclipsing — accre- 
tion, accretion disks — X-rays: binaries 



1. Introduction 

X-ray spectroscopy of low-mass X-ray binaries (LMXBs) provides quantitative constraints on 
the complex and poorly understood process of mass transfer through accretion disks. The discrete 
spectral features that are produced as the central continuum emission propagates through the 
circumsource material are sensitive to the physical conditions and geometry of the accreting plasma. 
High-resolution spectroscopy has been shown to be a powerful tool in empirically constraining the 
accretion process in LMXBs; analysis of the spectra obtained with the Chandra and XMM-Newton 
observatories have provided quantitiative measurements of the plasma parameters and accretion 
geometries in Cir X-l (Brandt & Schulz 2000) and EXO 0748-67 (Cottam et al. 2001). 

In this letter we present a preliminary analysis of the 4U 1822—37 spectra obtained with the 
Chandra High Energy Transmission Grating Spectrometer (HETGS). The smoothly-modulated, 



nearly color-independent partial eclipse {P = 5.57 hr) and the unusually hard X-ray spectrum of 
this source led White et al. (1981) and White & Holt (1982) to conclude that the central source is 
surrounded and partially obscurred by an accretion disk corona, formed from material evaporated 
off of an X-ray illuminated disk. Extensive modeling of the lightcurves by White & Holt (1982), 
Mason & Cordova (1982b), and Hellier & Mason (1989) constrained the geometric parameters of 
the system (i = 75 — 85°, .Rdisk = 5 — 7.3 x 10 10 cm, and -Rcorona ~ 2.5 — 4 x 10 10 cm), and suggested 
the presence of a phase-dependent vertical structure along the edge of the disk. The X-ray spectra 
obtained with previous low-resolution instruments have been extremely difficult to model; the data 
require a complex continuum and a broad Fe K line, however broad residuals at soft X-ray energies 
suggest the presence of additional unresolved emission or absorption structure (White et al. 1981; 
White, Kallman, & Angelini 1997; Parmar et al. 2000). 

With the high spectral resolution of the HETGS we are able to resolve discrete, narrow line 
emission from a range of abundant ions. As we describe below, detailed analysis of these features 
and their orbital phase variations suggest the presence of photoionized material on the inside edge 
of a bulge located at the point of impact with the accretion stream, and fluorescing material in the 
inner disk illuminated by scattered emission from the extended central corona. 



2. Data Reduction 

4U 1822-37 was observed on 2000 August 23 with the Chandra HETGS (Canizares et al. 2000) 
in the standard ACIS-S-HETG configuration. The observation began at 16:20:37 UT and lasted for 
approximately 40 ks covering two full binary periods. The source is bright with a phase-averaged 
flux of 4.8 x 10 -10 ergs/s/cm 2 in the 0.7 — lOkeV band. The zero order image is therefore severely 
piled-up, making it unusable for spectroscopy. However, in the combined positive and negative 
first order MEG and HEG spectra the phase-averaged count rates are only 2.8ct/s and 2.0ct/s 
respectively, and the spectra show no evidence of pile-up problems. 

We used the Chandra Interactive Analysis of Observations software (CIAO) to extract the 
spectra and generate the aspect-corrected effective area curves. However, due to limitations with 
the CIAO 1.1 version that was available the time, we were unable to generate the necessary files 
for time intervals within the observation. Instead we extracted the spectra for specific phases 
using custom-built software that was calibrated with the publicly available Capella data (see Be- 
har, Cottam, &: Kahn 2001). In order to utilize the latest effective area models, all global fitting 
was performed using the CIAO products within XSPEC (version 11.0: Arnaud 1996). The phase- 
dependent analyses were performed using our own software. For both procedures the positive and 
negative first order spectra were combined to maximize the statistical quality of the data. 



3. Observational Features 

3.1. Discrete Emission 

In Figure 1 we show the first order MEG and HEG spectra for 4U 1822—37 integrated over 
the total observation. Discrete emission lines are clearly visible from hydrogen-like and helium-like 
ions for a range of abundant elements including the Lya lines of Fe XXVI , Si XIV , Mg XII , Ne X , 
and O VIII , and the He-like complexes of Mg XI , Ne IX , and O VII . There is also a bright Fe Ka 
fluorescence line. A narrow radiative recombination continuum (RRC) from Ne IX is detected 
indicating that photoionization followed by recombination is the dominant excitation mechanism 
(Liedahl & Paerels 1996). We find no obvious absorption features in the spectra. 



3.2. Phase Dependence 

The observed spectral features are highly phase dependent. We divided the folded lightcurve 
(see Figure 2) into four phase bins centered on (ft = 0.00, 0.25, 0.50, and 0.75 and collected the 
spectra for each bin. In Figure 3 we show the MEG spectrum in the wavelength range of 5 to 15 A 
for the four phase bins. In Figure 4 we show the HEG spectrum in the wavelength range of 1.3 to 
2.5 A. The Ne X, Ne IX, Mg XII, Mg XI, Si XIV, and Fe XXVI emission lines are bright during 
phase (ft = 0.25. During phase (ft = 0.50 these lines are visible, but considerably weaker. A weak 
Mg XII emission line is visible during phase (ft = 0.0, and a weak Si XIV line is visible during phase 
(ft = 0.75. The oxygen region in not shown here since the phase-binned count rates are too low for 
statistically significant analysis. The phase-dependence of the recombination spectra suggests that 
the photoionized material comes from a localized region on the disk. 

The recombination lines show a small phase-dependent Doppler shift. It is difficult to determine 
precisely the centroid and width of the line emission, given the poor statistics in the phase-binned 
spectra. However, during phase (ft = 0.25 the emission lines are consistent with ^Doppler = (0 ± 
150) km/s. During phase (ft = 0.50 the emission lines show a slight blue shift with v Doppler = 
(280 ± 150) km/s. We see no velocity broadening of the lines with an upper-limit of At; < 300 km/s. 

The Fe Ka line has a very different phase-dependence. The line is bright during phases 
(ft = 0.0, 0.25, and 0.50. During phase (ft = 0.75 the line flux decreases noticably. At all phases the 
line centroid is consistent with 1.935 A, corresponding to transitions in cold, near-neutral material. 
We measure an upper limit to the line broadening of Av < 1500 km/s. 



3.3. Helium-like Series 

The He-like series of O VII , Ne IX , and Mg XI exhibit bright intercombination lines, weak 
resonance lines and essentially no forbidden lines (see Figure 5 and the insert on Figure 1). The 



ratio of the intercombination and forbidden emission to the resonance emission, G = (i + f)/r, is 
consistent with a purely recombining plasma; For O VII and Ne IX , where the lines are cleanly 
resolved, the measured ratio is G = 4.6 ±1.8 for O VII and 4.7 ± 1.5 for Ne IX , while the predicted 
ratios for an optically thin photoionized plasma are 4.0 and 3.7 respectively (calculated using 
HULLAC; Bar-Shalom et al. (1998)). The Mg XI series lies in the middle of a complex instrumental 
absorption feature making it difficult to quantify the line ratios. An enhanced intercombination line 
is often indicative of a high-density plasma. However, in 4U 1822—37, where there is a significant 
UV flux (3 x 10 -11 ergs/s/cm 2 from 1200 to 6400 A: Mason & Cordova 1982a), photoexcitation from 
the 2 3 S level to the 2 3 P level competes with the spontaneous decay of the 2 3 S" line to the ground 
state. Using the UV flux reported in Mason & Cordova (1982a) we calculated the photoexcitation 
rate according to the procedure described in Kahn et al. (2001). We calculated a dilution factor of 
~ 0.8 for the source flux based on the relative geometry of the corona and the accretion disk. For 
O VII and Ne IX we find the UV photoexcitation rate to be a factor of ~ 4 x 10 6 and ~4x 10 5 
larger than the rate of radiative decay. The ratios in the He-like series are therefore insensitive to 
the plasma density. 



4. Discussion 

These spectroscopic features can be used to locate and describe the line emitting regions in 
4U 1822—37. The phase-dependence of the recombination line intensity and the lack of absorption 
features are consistent with emission from an X-ray illuminated, optically thick region on the disk 
localized between binary phases (ft = 0.75 and (ft = 1.00. During phase (ft = 0.25 and part of (ft = 0.50, 
we see emission lines from the X-ray illuminated inside edge of the material. During phases (ft = 0.75 
and (ft = 0.00 this emission region is essentially blocked from view. This is illustrated in Figure 6. 
The inferred geometry is insensitive to our choice of phase bins. The observed velocity structure in 
the recombination lines is consistent with this scenario; emission lines from material localized near 
the edge of the disk at phase (ft = 0.85 will show very little velocity broadening, but will appear 
Doppler-shifted by the vector sum of the orbital and Keplerian velocities. At phase (ft = 0.25 
the lines will appear blue-shifted by ^Doppier ~ 140km/s and at phase (ft = 0.50 blue-shifted by 
^Doppier ~360km/s. 

The phase location of the emission region is consistent with the point of impact of the accretion 
stream and suggests that we are seeing line emission from the bulge that is expected to form in 
the shock- heated, colliding material (Lubow & Shu 1976; Livio et al. 1986). If the material is 
optically-thick, the bulge is expected to extend above the disk and downstream along the edge of 
the disk as it adiabatically cools (Armitage & Livio 1998). Viewed at high inclination angles, this 
would produce the wide shallow dip observed in the HETGS lightcurves and in those from previous 
observations. 

We can estimate the angular size of the material from the flux in the line emission by assuming 
that the photons above the ionization edge that intercept the material are absorbed and ultimately 



reemitted as recombination line photons. Expressing this in terms of the number flux, as iV con ti nU um * 
Jl/(4vr) = Absorbed = ^rccombincd, and iV rccombincd * r] Vmc = N Vmc , where T] is the efficiency of cascade 
through the line we can calculate the solid angle subtended, f2. Using the spectrum from the 
4> = 0.25 phase bin we find a solid angle of 0/(4vr) = (1.8 ± 0.7) x 10" 2 . For the 4U 1822-37 
geometry this corresponds to an illuminated surface area of (1.0 ± 0.4) x 10 21 cm 2 , or a vertical 
scale height of ~ (3.3 ± 0.7) x 10 10 cm. 

The presence of a bulge at the stream impact point in 4U 1822—37 was first suggested by 
Mason et al. (1980) to explain the sinusoidal modulation in the optical lightcurve. Using detailed 
modeling of the X-ray, UV, optical, and IR lightcurves White &; Holt (1982), Mason & Cordova 
(1982b), and Hellier & Mason (1989) all inferred a phase-dependent vertical structure along the 
edge of the disk that is dominated by a bulge around phase <f> = 0.85, which extends from <j) = 0.70 
to <fi = 0.0 with a maximum vertical height between 0.6 and 1.6 x 10 10 cm depending on the model. 
This is consistent with the phase location and the vertical height that we measure here to within 
a factor of two. 

We can use the details of the HETGS spectra to further constrain the physical parameters of 
the emitting material. The width of the RRC feature is a direct measure of the plasma temperature. 
For the cleanly resolved Ne IX RRC we find a temperature of 13=L7eV = (1.5=L0.8) x 10 5 K. We can 
estimate the plasma density of the emitting material using the emission measure, EM = f dVn^. 
The emission measure is determined by dividing the line luminosity by the line powers, which are 
calculated using HULL AC (Bar-Shalom et al. 1998). For all ions except oxygen, we measure EM ~ 
5 x 10 56 (D/2.5kpc) 2 cm" 3 . For O VIII and O VII we measure EM ~ 2.5 x 10 55 (D/2.5kpc) 2 cm" 3 , 
which could imply a relative underabundance of oxygen by a factor of ~ 20. Estimating the volume 
by taking the vertical height that we infer here, the azimuthal extent from previous modeling of the 
lightcurves, and assuming that the emission measure is dominated by material within a Thompson 
depth of r = 1, we find an average electron density of ~ 1 x 10 11 cm" 3 . Such a high density would 
have been discernable in the He-like line ratios if the UV photoexcitation had not saturated the 
transitions. 

The persistence of the Fe Ka fluorescence emission throughout the binary phases with a de- 
crease in the intensity only during phase <p = 0.75 suggests that the fluorescing material comes from 
an extended inner region of the system, which is visible at all times, but partially blocked from 
view at that point by the bulge. The width of the line constrains the radial location of the emitting 
material. From the measured upper-limit to the velocity broadening we find a lower limit for the 
inner radius of R > 1 X 10 10 cm. According to models for the formation of the corona, it should be 
highly ionized and optically thick (White & Holt 1982). The observed Fe Ka line cannot originate 
in such a corona or else it would have a higher energy centroid and would appear broadened by 
Compton scattering in the hot (~ 1.8 keV) coronal material. 

We can estimate the relative geometry of the fluorescing material and the illuminating central 
source using the flux in the fluorescence emission line. Calculating the illuminated solid angle 



using the same assumptions as before, but now with the fluorescence yield of neutral iron, we find 
J7/(27r) = 0.11 ± 0.3. This is much larger than expected for an embedded point source illuminating 
a disk, unless the edge of the disk flares to scale heights comparable to the bulge height. In order to 
illuminate such a large solid angle, the central source must be extended. However, the fluorescing 
material must be even more extended; an extended source laying over an illuminated region of 
comparable size produces a solid angle of 0, ~ 2tt, which is much larger than we observe here. 
The cold clumps formed by instabilities in an accretion stream that skims over the disk surface 
(Frank, King, & Lasota 1987) might provide sufficient solid angle. Alternatively, the disk around 
an embedded spherical central source subtends a solid angle of Q/(2tt) < 0.25 (Chen & Halpern 
1989). A covering fraction of 0.11 occurs when the inner radius of the illuminated material is about 
twice as large as the radius of the spherical source. In 4U 1822—37 this would give a lower-limit 
to the radius of the corona of i? C orona > 5 x 10 9 cm, which is only a factor of ~ 5 smaller than the 
radius determined by modeling the lightcurve. 

Support for the work at Columbia University was provided by the NASA through Chandra 
Award number SAO-G00-1099X issued by the Chandra X-ray Observatory Center under contract 
NAS8-39073. Work at LLNL was performed under the auspices of the U.S. Department of Energy 
by the University of California Lawrence Livermore National Laboratory under contract No. W- 
7405-Eng-48. 
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Fig. 1. — HETGS spectra of 4U 1822—37. The thin and thick histograms correspond to the MEG 
(m = ±1) and HEG (m = ±1) spectra, respectively. The solid lines represent a simple fit to the 
continuum emission and are intended to indicate the instrumental features in the effective area. 
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Fig. 2. — Lightcurve for the HETGS observation in A0 = 0.01 phase bins. The thin and thick 
histograms show the MEG (m = ±1) and HEG (m = ±1) lightcurves, respectively. 
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Fig. 3. — MEG (m = ±1) spectrum of 4U 1822—37 for four phase bins centered on 
0.0,0.25,0.50,0.75. 
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Fig. 4. — HEG (m = ±1) spectrum of 4U 1822—37 for the four phase bins. 
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Fig. 5. — Ne IX complex in the MEG (m = ±1) spectrum. The resonance (R), intercombination 
(I), and forbidden (F) lines are labeled. 
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Fig. 6. — Illustration of the inferred geometry in the 4U 1822—37 system. This illustrates the 
phase-dependent obscuration of the central source by the bulge, as well as the phase- varying aspect 
of the irradiated bulge surface. 



